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Reply to comment on "Competing Interactions, 
the Renormalization Group and the Isotropic- 
Nematic Phase Transition" 

The focus of our work [l| was to identify conditions for 
the presence of an isotropic-nematic phase transition in 
the context of a generic system with isotropic compet- 
ing interactions. By taking into account non trivial an- 
gular momentum contributions from the interaction, we 
found a second order isotropic-nematic phase transition 
at mean field level, which becomes a Kosterlitz-Thouless 
one Q when fluctuations are taken into account. 

In his comment, Levin criticizes our results by showing 
that the low temperature fluctuations of a stripe phase 
in 2d diverge linearly in the thermodynamic limit. His 
analysis is restricted to the stripe phase and, contrary 
to what is suggested in the comment, does not apply to 
the central result of our letter which is the existence of 
an isotropic-nematic phase transition. In fact, as clearly 
anticipated by us in the letter [l[ , the corresponding anal- 
ysis of the fluctuations of the nematic order parameter 
displays a logarithmic divergence leading to a low tem- 
perature phase with quasi long range order. 

In our model, despite of the involved calculations, it 
is straightforward to understand this fact. Introduc- 
ing the nematic order parameter Qij = a (fiihj — 
(where hi = (cos 8, sin 9) is the director field) through a 
Hubbard-Stratonovich transformation, it is possible to 
decouple the quartic <p terms. Integrating out the <j) 
field, we obtain the following long wavelength effective 
free energy for the nematic order parameter: F(Q) = 
a 2 /2 Tr{Q 2 ) + o 4 /4 Tr(Q 4 ) + p/4 Tr(QDQ) + . . ., where 
the symmetric derivative tensor Dij = V^Vj and 04 
and p are temperature dependent coefficients given in 
terms of the parameters of the original model. At mean 
field, the last term is zero, and we find a = y— 02/04 for 
02 < 0, going continuously to a = for 02 > 0. Note that 
any global rotation of the order parameter costs no en- 
ergy. Therefore, parametrizing the order parameter by a 
modulus and an angle, the long wavelength angle fluctu- 
ations 6{x) dominate the low energy physics. Computing 
the free energy at lowest order in the derivatives of the an- 
gle fluctuations we find AF = pa 2 J d 2 x\V9\ 2 , where AF 
is the excess of free energy relative to the saddle point 
value. Therefore, the free energy of fluctuations corre- 
sponds to that of the XY model. The only difference with 
the usual vector orientational order is that the system 
should have the symmetry 9 — > 9 + n modifying the vor- 
ticity of the topological defects. Thus, one finds for the 
angle fluctuations < 9(x)9(x') >~ ln[fco(ai ~ x')\, which 


in turn lead to an algebraic decay of the order parame- 
ter correlations. In an extended paper we will show the 
explicit dependence of the Frank constant K(T) = pa 2 
with the parameters of our model fcoj m > uq and Ui. 
The conclusion is that the isotropic-nematic transition 
takes place by the Kosterlitz-Thouless mechanism of vor- 
tex (disclination) unbinding This result agrees with 
the predictions of the well known KTHNY theory, that 
predicts the same kind of transition in a layered two- 
dimensional system |3| . The main difference between our 
work and the KTHNY theory is that while KTHNY begin 
the analysis from an elastic energy valid at low temper- 
atures already in the crystal phase, and explicitly add 
a term to take into account topological defects, we ap- 
proach the transition form the disordered phase, allowing 
the possible emergence of spontaneous symmetry break- 
ing. In this way, we make contact with more microscopic 
parameters, clarifying in some way the role of competing 
interactions in the development of orientational phases, 
contrary to what is suggested in the comment. It is well 
known that Ginzburg-Landau functionals similar to the 
one explored by us can be obtained from more micro- 
scopic interactions like those in ultrathin magnetic films 
[J] and copolymers in 2d [a]. 

In conclusion, we have shown that the model intro- 
duced in equations (I) and (7) of Ref. [l[ undergoes 
an isotropic-nematic phase transition in the Kosterlitz- 
Thouless universality class. 
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